Observation of strong spin— orbital entanglement in Sr2Ru04 by spin-resolved ARPES 
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After a flurry of experimental activity [1-5], Sr2Ru04 has become a hallmark candidate for spin- 
triplet chiral p-wave superconductivity, the electronic analogue of superfluid "^He [6]. However, de- 
spite the apparent existence of such a pairing, some later experiments [7-9] do not fully support this 
conclusion, as they cannot be explained within a theoretical model using spin-triplet superconduc- 
tivity alone [10]. A resolution might come from the inclusion of spin-orbit coupling, which has been 
conjectured to play a key role in the low energy electronic structure [11]. Using circularly polarized 
light combined with spin- and angle-resolved photoemission spectroscopy, we directly measure the 
value of the effective spin-orbit coupling to be 130=b30meV. This is even larger than theoretically 
predicted and comparable to the energy splitting of the dxy and dxz,yz orbitals around the Fermi 
surface, resulting in a strongly momentum-dependent entanglement of spin and orbital character. 
As a consequence, the classification of the Cooper pairs in terms of singlets or triplets fundamentally 
breaks down, necessitating a description of the unconventional superconducting state of Sr2Ru04 
beyond these pure spin eigenstates. 



With the discovery of topological insulators [12], as 
well as relativistic Mott-insulating behavior in 5d transi- 
tion metal oxides [13, 14], there has been a renaissance 
in the study of spin-orbit (SO) coupling and the novel 
phenomena it might induce in complex materials. By 
mixing the canonical spin eigenstates, the relativistic SO 
interaction might also offer a resolution to conflicting 
experiments regarding the nature of the superconduct- 
ing pairing in Sr2Ru04. In fact, SO coupling's major 
role in determining the electronic band structure - as 
indicated by ab initio density functional theory results 
[11] - suggests that its consideration may indeed be im- 
portant when describing superconductivity as well; this 
could play a fundamental role beyond simply lifting the 
degeneracy of competing pairing states [11, 15-18]. 

However, thus far, the experimental investigation of 
SO coupling's effects on the electronic structure of 
Sr2Ru04 has been limited to the comparison of band 
calculations against angle-resolved photoemission spec- 
troscopy (ARPES) [11, 19, 20] - no success has been ob- 
tained in observing experimentally either the strength 
of SO coupling or its implications for the mixing be- 
tween spin and orbital descriptions. Here we probe 
this directly by performing spin-resolved ARPES with 
circularly-polarized light, by using the angular momen- 
tum inherent in each photon, along with electric-dipole 
selection rules, to generate spin-polarized photoemission 
from the SO mixed states. Combined with a novel spin- 
and orbitally-resolved ab initio-based tight-binding (TB) 
modelling of the electronic structure (Methods), these 
results demonstrate the presence of a non-trivial spin- 



orbital entanglement over much of the normal-state Fermi 
surface, i.e. with no simple way of factoring the band 
states into the spatial and spin sectors. 

In Sr2Ru04 the calculated effective SO coupling is 
comparatively small (Ceif ^ 90 meV at the F point) with 
respect to the bandwidth (~ 3eV) of the Ku-t2g orbitals, 
which define the a, /3, and 7 conduction bands. Never- 
theless its influence always becomes important whenever 
bands would be degenerate in the absence of SO, either 
by symmetry or accidentally. This happens at several 
places in the three-dimensional (3D) Brillouin zone, as 
demonstrated in Fig. l(a,b) where we show a compari- 
son of the ab initio-TB band structure and Fermi surface 
calculated both with (color) and without (black) SO cou- 
pling included (Methods). In the absence of SO, by sym- 
metry the dxz and dy^ bands would be degenerate along 
the entire kz momentum path from F to Z [Fig. 1(a)]. 
Additionally, there are accidental degeneracies along the 
kz = path from F to X, where the bands of dxz,yz and 
dxy character all cross at momenta near (27r/3, 27r/3) - 
the exact location of which varies with k^ but often occurs 
at the Fermi level [Fig. l(a,b)]. At all these locations SO 
coupling naturally leads to a splitting [Fig. l(a,b)] and 
mixing of the orbital character [Fig. 1(c)] for all three 
bands. 

Interestingly, the effects of SO coupling are not limited 
to the regions around the non-relativistic degeneracies as, 
despite the large bandwidth, the Ru-t2g bands are often 
separated by energies small compared to the SO interac- 
tion. The predicted importance of the SO interaction can 
be directly visualized via the expectation value of /-s from 
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FIG. 1: Effect of SO coupling on the band dispersion and Fermi surface of Sr2Ru04. (a) Band structure along the 
high-symmetry directions, and (b) kz = Fermi surface, calculated both without (thin black) and with (thick, color-coded 
to show (/ • s)) the inclusion of SO coupling; at the F point, the latter gives rise to a (ejj 90meV splitting [note that 
Z= (0,0,7r/c), F= (0,0,0), M= (7r/a,0,0), X= (vr/a, vr/a, 0)]. (c,d) 3D Fermi surface sheets color-coded to show (c) orbital 
character and (d) the expectation value {I ■ s); the extended zone scheme is used to illustrate the conventional Brillouin zone in 
the body-centered tetragonal unit cell of Sr2Ru04. The energy and momentum location of the spin-resolved ARPES spectra 
shown in Fig. S3 is marked in yellow in panels (a) and (b). 



our ab initio-TB modelling, with / and 5 being the orbital 
and spin angular momentum operators. A non-zero value 
of (l-s) indicates complex orbital eigenstates that can be 
entangled with the spin. In this case, the wavefunction 
cannot be factorized into independent spin and orbital 
parts, as would be possible for a fully quenched angular 
momentum (for which (l-s) =0). The calculated (l-s) 
is shown in Fig. 1 for the high-symmetry dispersion (a), 
kz = Fermi surface (b), and around the 3D Fermi sheets 
(d). This suggests SO coupling is important in Sr2Ru04 
on almost the entire three-sheet Fermi surface [21]. 

In order to probe the resulting internal spin-orbital 
structure of the band states, we turn to spin-resolved 
ARPES with circularly polarized light: with this tech- 
nique the circular polarization of the light couples to the 
angular momentum of the states measured at a given 
k point, while the spin is resolved directly. A similar 
approach, albeit without the angular and energy resolu- 
tion needed to resolve states near the Fermi level, has 
been used previously to generate spin-polarized photoe- 
mission from materials without a net magnetization, such 
as GaAs [22] and Ca2Ru04 [23]. Here, exploiting the 
electron-dipole selection rules for photoemission from ini- 
tial states carefully selected via spin-resolved ARPES, 
we will - for the first time - directly probe the internal 
SO structure of the wavefunction for these critical near- 
Fermi-level states. This study will be done at ~40K, 
thus in the normal state well above Tc — 1.5K. 

To apply this technique on Sr2Ru04 we study the SO 
splitting at the F point, k = (0, 0, 0); this location is high- 



lighted in Fig. l(a,b). As explained below, this selects 
both the experimental geometry and initial-state wave- 
functions that are the most straightforward to under- 
stand, facilitating the direct measurement of both the SO 
interaction strength and the complex nature of the wave- 
function. At this /c-point, non-relativistic band structure 
calculations predict two degenerate bands of dxz and dyz 
character, with the dxy band far enough away that it 
can be ignored (i.e., at about 1.8-2.3 eV higher binding 
energy, depending on the kz value). Here SO breaks the 
degeneracy by hybridizing these bands to form two states 
with a splitting of (eff ^ 90 meV: a lower binding-energy 
state with z-components of orbital and spin angular 



momentum parallel 



and a higher binding- 



energy state where they are antiparallel 
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Optical selection rules for the initial-to-final-state exci- 
tation with circularly polarized light dictate that both 
= zbl and Am^ = ±1. For d orbit als the change in 
^ will strongly favour the d ^ p over d ^ f transitions, 
owing to the cross-section at the photon energies used 
(24 and 56 eV); the change in mi will depend on the cir- 
cular polarization of the photon being right (0) or left 
(G). When a (0) photon is absorbed by the lower 

binding-energy parallel d'^^^ , d^^^ ^ state, mi^ must in- 
crease (decrease) by one; but since \mi^ \ = 2 is forbidden 
in the strongly favoured p transition, electrons from the 
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FIG. 2: Spin-resolved ARPES results on Sr2Ru04 at the F point, (a) Measured polarization asymmetry of the 
photoemitted electrons, and (b-d) corresponding spin-resolved ARPES intensities relative to the x, and z crystal axes 
obtained with right (©) or left (0) circularly polarized 24 eV photons at the F point, as highlighted in Fig. 1. Note that, for 
the z direction (d), the intensity from each underlying state is shown, corrected for light incident at 45° with respect to the 
spin-orbit quantization axis, as detailed in the Supplementary Information. Error bars represent statistical uncertainty based 
on number of counts in the Mott polarimeters, plotted at 95% confidence, shown together with locally- weighted scatterplot 
smoothing fits (Supplementary Information). 



d^-i^ (^+1^) ^^^^ degenerate state will dominate, 

resulting in an effective (t^) spin polarization. Simi- 
larly, photoemission from the higher binding-energy an- 
tiparallel state (it.\^5^+i^^ using (0) light will result 
in photoemission with the opposite spin polarization, 

When observed with regular spin-integrated ARPES, 
the energy distribution curve from the F point shows 
both predicted states as a single broad feature with width 
~ 400 meV (Supplementary Information); however, by 
using circularly polarized light and observing the spin- 
polarization of the photoelectrons, it is possible to dis- 
tinguish them. The experiment is repeated twice, with 
both helicities of light, and the results combined to cal- 
culate the photoelectron polarization asymmetry which 
eliminates possible experimental artefacts (Supplemen- 
tary Information) . This polarization asymmetry is shown 
in Fig. S3(a): it is zero for the x and y directions, and 
shows a clear wiggle as a function of energy in the z di- 
rection. The non-zero values indicate that the photoelec- 
trons have a photon-helicity-dependent spin-polarization 
only in the z-direction, while the binding-energy depen- 
dence is due to the origin of the photoelectrons shifting 
from parallel to antiparallel states. By plotting the inten- 
sities corresponding to the observed photoelectron polar- 
ization asymmetry for each spatial dimension. Fig. S3(b- 
d), we can directly resolve these states. In Fig. S3(d), the 
z direction, they become visible as two energy-split fea- 



tures: 



is the state that photoemits with 
light and with light, i.e. the feature detected 



in la 



whereas 



is detected in I^a^^q^. 

Along the x and y directions in Fig. S3(b,c), however, 
the spectra match the spin-integrated intensity since the 
photoelectrons from both states have {sx) = {sy) =0 
for both helicities of light. The splitting in the z direc- 
tion is observed with both 24 eV and 56 eV photons and 



its magnitude is 130 ± 30 meV (Supplementary Informa- 
tion) , showing good agreement or a possible enhancement 
compared to the predicted SO splitting of C^eff ~90meV. 
Most importantly, the existence of these two states, from 
which spin-polarized photoemission can be generated us- 
ing circularly polarized light in the z direction only, is 
clear experimental evidence of the importance of SO cou- 
pling in Sr2Ru04 and of its consequences for the complex 
nature of the normal-state wavefunctions. 

As discussed below, the most important of these con- 
sequences is the strong, momentum-dependent, spin- 
orbital entanglement of the eigenstates around the Fermi 
surface. This is illustrated in Fig. 3 by plotting the pro- 
jection of the Bloch wavefunctions at the Fermi energy 
onto the Ru-d orbital basis at different momenta (Meth- 
ods). The resulting projections are color-coded by the ex- 
pectation value of the spin operator {sz)(e,(/)) fo^ half 
of the Kramers-degenerate pair (blue='f, red=|). Along 
the edges of the Brillouin zone (M— X) where the bands 
are well-separated, we find particularly in the a band (far 
right panel in Fig. 3) that the orbitals do not show strong 
entanglement: each orbital projection is associated with 
a single expectation value (color) of the spin operator. 
We also notice that in these areas the /3 and a bands are 
of pure dxz,yz orbital character, and the 7 band of dxy 
(Fig. 3). At these locations the wavefunction could be 
well approximated by the usual description as a product 
of independent spatial and spin components: 



^(k,<7)=^(k) 



(1) 



where (fCk) and are the spin and orbital eigenstates, 
and a the spin index. However, close to the zone diago- 
nal, e.g. near the intersections of the Fermi sheets with 
F— X, this is not the case. Here we find strong orbital mix- 
ing for all bands and, especially in the 7 and /3 bands, 
also strong entanglement between orbital and spin char- 
acter: the orbitals are no longer associated with a uni- 
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FIG. 3: Momentum dependent spin— orbital entanglement of the Sr2Ru04 eigenstates around the Fermi surface. 

Momentum- dependent Ku-d orbital projection of the wavefunction for the /3, 7, and a bands at selected momentum locations on 
the 3D Fermi surface. The surface color represents the momentum-dependent Sz expectation value along the direction defined 
by the spherical {0,(j)) angles, {sz)(e,(j)) (see Methods for definition); as indicated by the color scale at upper left, blue/red 
correspond to spin t/i for one state of the Kramers-degenerate pair (with the other state - not shown - having opposite spin 
polarization, consistent with the absence of net spin polarization at any k point in Sr2Ru04). The strongly mixed colors on 
some of the orbital projection surfaces indicate strong, momentum-dependent spin-orbital entanglement. 



form spin value (color) and, on the contrary, the latter 
can vary from fully up to fully down along a single orbital 
projection surface. Meaning that, at these locations, the 
wavefunction cannot be expressed as the product of in- 
dependent spatial and spin factors; instead, we need the 
more general expression: 

^(k, a) = ctn(k)0r" + c;^;(k)0f " , (2) 

with a being the pseudo-spin index, and c^^^ the prefac- 
tors of the momentum-dependent spin-orbital-entangled 
eigenstates. Eq. 2 further illustrates the nature of the 
SO-induced entanglement: flipping the spin forces also a 
change of the orbital character. We note that, due to the 
nature of the band structure in Sr2Ru04, this entangle- 
ment is strongly momentum dependent in both k\\ and 
kz, despite the extremely weak kz dispersion along the 
Fermi surface. 

A similar momentum and orbital dependence of the 
spin expectation value is responsible - in topological 
insulators - for the complex spin-texture of the Dirac 
fermions [12, 24]. In Sr2Ru04, beyond the normal-state 
properties, it directly affects the description of supercon- 
ductivity, as revealed by the inspection of the Cooper 
pair basic structure. Cooper assumed the two-particle 
wavefunction describing a superconducting electron pair 
to be of the form ^(ri, ai, r2, (72) = (p{ri - r2 ) ' 
with zero total momentum and the spin part being ei- 
ther singlet (total spin 6' = 0) or triplet {S = l) [25]. This 
allows one to classify superconductors as a realization of 
singlet or triplet paired states. However, a fundamental 
assumption of this description is that one can write the 
wavefunction of each electron as a simple product of spa- 
tial and spin parts, which is not possible in the case of 



strong mixing between (p^Ck) and ^iCk). Additionally, 
because of the strong 3D /c-dependence of this entangle- 
ment in Sr2Ru04, any transform to pseudo-spin would 
also necessarily be /c-dependent - negating the possibil- 
ity of using the regular description under a pseudo-spin 
basis as might be done, e.g., when describing supercon- 
ductivity in the heavy-fermion Ce compounds [26, 27]. 
As a consequence, the classification of Cooper pairs in 
terms of singlets or triplets fundamentally breaks down 
for Sr2Ru04. This is shown in Fig. S5 of the Supple- 
mentary Information, which presents the spin-eigenstates 
available to a pair of electrons with zero total momentum 
as obtained from the expectation value ( si • S-k )• While 
the familiar singlet and triplet states are seen off the zone 
diagonal for the a band (with (sk-S-k) = —3/4 and 1/4, 
respectively), they are not available for either the /3 or 7 
bands, whose spin-eigensystem consists of a doublet and 
two singlets or - depending on the location in /c-space - 
two doublets. 

Our findings mark a deviation from a pure spin-triplet 
pairing for Sr2Ru04, since the only portion of the Fermi 
surface that might support it is contained within the 
a pocket. This could explain a number of experi- 
mental observations at variance with a spin-triplet sce- 
nario, such as the extreme sensitivity to field angle of 
both the magnetic-field-induced second superconducting 
phase transition [7] and also the suppression of the ab- 
plane upper critical field [8]. These provide evidence for 
an additional magnetic anisotropy in the superconduct- 
ing state, of which SO coupling would be the natural 
source. Our study suggests that the superconductivity 
in Sr2Ru04 is yet more unconventional than has been 
assumed so far: provided that the paring interaction is 
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relatively weak, the entanglement of spin and orbit at the 
single-particle level becomes an essential feature to take 
into account when describing the superconducting pair- 
ing. In this regard, it would be interesting to verify what 
of the chiral p-wave superconductor phenomenology [1- 
5], and apparent conflict in experimental evidence [7-9], 
would remain when re-evaluated in terms of entangled 
single-particle eigenstates. 

METHODS 

Spin-resolved ARPES experiments. The mea- 
surements were performed on the normal state of 
Sr2Ru04 at the Swiss Light Source using the COPHEE 
endstation on the Surface and Interface Spectroscopy 
beamline, with ^lOOmeV energy resolution and ~1° an- 
gular resolution (the latter is equivalent to ^O.lir/a 
at 24 eV and ^0.2 7r/a at 56 eV photon energy), which 
allows the selection of electrons in energy-momentum 
space. These electrons are deflected into an array of Mott 
detectors, where the spin polarization can be resolved 
in all three spatial dimensions. High-quality Sr2Ru04 
single crystals with no 3K phase contamination [28] 
were cleaved in situ at 40 K and base pressures in the 
10~^ mbar range. Owing to the relatively high tempera- 
ture and pressure, no features associated with the recon- 
structed surface were detected [28]. All measurements 
were repeated using the beamline's elliptically-polarizing 
undulator to generate both right (0) and left (o) circu- 
larly polarized light, and the photoelectron polarization 
asymmetry was calculated using the geometric for- 
mula [29] to eliminate the effects of circular dichroism 
(Supplementary Information): 

^*^^Motts = — / / • (3) 

Here 5'Motts is the empirically determined measure of the 
efficiency of a Mott spin-detector pair (the so-called Sher- 
man function), and / is the intensity for the up {U) and 
down (D) halves of the Mott pair measured with circular 
plus (0) and minus (0) light. This formula was used for 
each of the Mott polarimeter pairs and the resulting po- 
larization asymmetry was then translated into sample co- 
ordinates according to the experimental geometry. With 
the effects of dichroism eliminated, this photoelectron po- 
larization asymmetry is equivalent to the true energy- 
and momentum-resolved photon-helicity-dependent elec- 
tron spin polarization. 

Electronic structure calculations. To theoret- 
ically study the effects of the SO interaction on the 
low-energy electronic states of Sr2Ru04, and in particu- 
lar the momentum-dependent spin-orbital entanglement, 
the bulk band structure was calculated using the linear 
muffin-tin orbital (LMTO) code, which was tested 
against that calculated with full-potential linearized 



augmented-plane-wave method using WIEN2k. The 
parameters were extracted from the LMTO results 
using the order-N muffin-tin orbital method [30] on the 
minimal basis of oxygen 2p and ruthenium Ad orbit als. 
Atomic SO coupling was then added as a local term in 
the Hamiltonian, the results of which show excellent 
agreement with the Fermi surface and band dispersion 
obtained from relativistic (i.e., SO included) calculations 
in WIEN2k. The orbitals in Fig. 3 depict the (l9, (/)) 
angular dependence of the Ku-d orbital projection of 
the Bloch wavefunctions for a given momentum at the 
Fermi energy. Shown is the surface defined by the 
equation r(l9, ^) = Er,r' ^r(6>, ^)^r' (6>, ^)(a[^^;,a.^^,^^), 
where i and r are the band and Ru-d orbital indexes, k 
is the momentum of the Bloch eigenstate, and Z are the 
cubic harmonics. The surface color corresponds to the 
momentum-dependent Sz expectation value along the di- 
rection defined by the spherical (6>, 0) angles, {sz)(^e,(j)) = 

with only one of the Kramers degenerate states shown. 
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SUPPLEMENTARY INFORMATION 
Components of {I ■ s) around the Fermi surface 

In order to investigate the coupling between the spin 
and orbital angular momentum near the Fermi energy 
in Sr2Ru04 we plot the expectation value of their dot 
product, ( / • s), for the states around the Fermi surface 
sheets. In a system without spin-orbit (SO) coupling, 
the independence of these two vectors would cause this 
expectation value to be zero everywhere. In Fig. Sl(b) 
we show ( / • s), contrasted against the orbital character 
in (a), as presented in the main text. In (a) we see how 
the /3 and a bands are of mostly dxz and dyz character, 
while the 7 band is of mostly dxy character; the excep- 
tion is near the anticrossing of the P and 7 bands close 
to (^, ^) where the strongest mixing occurs and their 
orbital character changes. In (b) we see that, despite 
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FIG. SI: The Fermi surface sheets of Sr2Ru04 colorized to 
show (a) orbital character, as well as (b) {l-s) and its compo- 
nents (c-e) along the three crystal axes (for the latter the color 
scale is shown at the bottom). The extended zone scheme is 
used to illustrate the conventional Brillouin zone in the fun- 
damentally body-centered tetragonal unit cell. 



the apparent mixing only at the anticrossing noted in 
(a) , spin and orbital angular momentum remain coupled 
around almost the entire Fermi surface, as evidenced by 
the non-zero value of (l-s). Interestingly, we find that 
(l-s) = only near the small portion of the 7 band 
which is of d^z/dyz character. However, in (c-e) we see 
that this is not because spin and orbit are independent 
at this location in momentum space; instead we see, in 
(e) for the z direction especially, that / and 5 are just as 
strongly coupled in this location but that (l-s) =0 only 
as a sum of non-zero components. This demonstrates 
that the spin and orbital angular momentum are coupled 
everywhere on the Fermi surface, for all three Ru-t2^ con- 
duction bands. 
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FIG. S2: ARPES on Sr2Ru04 at the F point; (a) shows an 
EDC integrated over ^0.03 ^ around the F point, as marked 
in the band map of panel (b). In (b) bands are labelled as 
they cross the Fermi level (with primes denoting reconstructed 
features [28]); the location of the band map is marked on the 
Fermi surface map shown in the inset. 



uum with pressure in the 10~^^ mbar range, using the 
21.22 eV He I spectral line, and resolutions of ~17meV 
and ~0.01 ^ angular equivalent. FigS2(a) shows an en- 
ergy distribution curve (EDC) taken at the exact same 
spot in 2D momentum space as that measured with spin- 
ARPES and presented in the main text, while FigS2(b) 
shows its location in momentum space on the measured 
band map and Fermi surface. This data demonstrates 
that the two states studied by spin- ARPES and presented 
in the main text have a fundamentally broad lineshape - 
even with a resolution approximately seven times smaller 
than the splitting their separation cannot be resolved 
without circularly polarized light and spin discrimina- 
tion. Note that this sample, cleaved at lower temperature 
and in better vacuum conditions, shows clear evidence of 
the reconstructed surface electronic structure [19, 28] (as 
seen by the folded a' band), while that presented in the 
main text did not. In either case no additional folded 
features will appear in the studied EDC at the F point, 
owing to the nature of the folding which maps the X 
point (which has no features in that energy range) to the 
F point - an additional benefit of performing the exper- 
iment at this location in momentum space. 



High resolution ARPES data 

In addition to the spin-resolved data presented in the 
main text, higher-resolution angle-resolved photoemis- 
sion spectroscopy (ARPES) experiments were performed 
at the University of British Columbia. In FigS2 we 
present some of this data taken at 5.2 K, under vac- 



Polarization asymmetry and the angle of incidence 
in spin- ARPES 

The geometric polarization asymmetry [29], defined as 
in the main text in terms of the intensities (/) from the 
up ([/) and down {D) detectors for circular right (0) and 
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FIG. S3: Spin-ARPES using 24eV circularly polarized light 
on Sr2Ru04 at the F point, (a-c) Spin-ARPES intensities rel- 
ative to each of the crystal axes obtained with right (0) or left 
(©) circularly polarized light; these are derived, according to 
Eq. S7, from the integrated spin data and measured polariza- 
tion asymmetry shown in panel (d) and (e). (f) Intensity from 
each underlying state, corrected for light incident at 45° with 
respect to the spin-orbit quantization axis and assuming no 
background, according to Eq. S9. Error bars represent statis- 
tical uncertainty based on number of counts in the Mott Po- 
larimeters, plotted at 95% confidence, with locally weighted 
scatterplot smoothing fits shown [31]. 



left (g) light, is given by the following expression: 



Motts 



j©j0 



t^tQ 



t®tO 



t^tQ 



(S4) 



which eliminates the effects of circular dichroism and any 
response imbalance that might exist between the two de- 
tectors in the Mott pair (the Sherman "function", ^Motts, 
is a measure of the efficiency of spin discrimination and 
is a constant, which we ignore here for simplicity). The 
effectiveness of Eq. S4 can be seen by modelling the mea- 
sured intensities as: 

1% = (PvDlm.ei 
1% = d^rol^uet 



(S5) 



Here the d and r terms account for the circular dichro- 
ism and detector response rates and /®^,e| = I^^ = Iqi 
(and its opposite, I^^^q^) are the fundamental photoe- 
mission intensities we are interested in, whose spin po- 
larization changes sign upon switching the helicity of the 
incident light. By substitution we see that the dichroism 
and detector efficiencies cancel out, and the polarization 



FIG. S4: Spin-ARPES using 56eV circularly polarized fight 
on Sr2Ru04 at the F point, and taken over a larger energy 
range than Fig. S3; this wider range shows that the polar- 
ization asymmetry goes to zero at the non-SO split higher 
binding-energy states, (a-c) Spin-ARPES intensities relative 
to each of the crystal axes obtained with right (0) or left 
(0) circularly polarized light; these are derived, according to 
Eq. S7, from the integrated spin data and measured polariza- 
tion asymmetry shown in panel (d) and (e). (f) Close-up in 
energy of the intensity from each underlying state, corrected 
for light incident at 45° with respect to the spin-orbit quanti- 
zation axis and assuming no background, according to Eq. S9; 
this panel is taken from an additional data set with narrower 
energy binning. Error bars represent statistical uncertainty 
based on number of counts in the Mott Polarimeters, plotted 
at 95% confidence, with locally weighted scatterplot smooth- 
ing fits shown [31]. 



asymmetry reverts to the simple form: 



(S6) 



The total intensity {It) can thus be combined with the 
measured polarization asymmetry to reconstruct the spin 
intensities for each spatial direction according to: 



W; = ^t(1 + P®)/2 
/e4,et = ^T(l-P®)/2. 



(S7) 



These are the spin intensities shown in Figs. S3 (a-c) and 
S4(a-c) for the two measured photon energies of 24 and 
56 eV, as constructed directly from the integrated spin 
(d) and polarization asymmetry (e). From these datasets 
two distinct features, split in the z direction only, are 
already apparent (we also note that, while similar in all 
respects, the wider binding energy range 56 eV data show 
that the polarization asymmetry vanishes at the non-SO 
split, high-energy oxygen states). However, owing to the 
experimental geometry (light incident at 45° relative to 
the analyzer), photons cannot be incident directly along 
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the z direction (additionally there is a node in the final 
state preventing the experiment from being performed in 
the second Brillouin zone with normal incident light and 
45° emission). The result is that with equal probabil- 
ity the photons will interact either with mi^ , resulting in 
spin-polarized photoemission from the states with paral- 
lel and antiparallel spin and orbital angular momentum 
(as described in the main text), or with m^^, resulting 
in non-polarized photoemission from both states (which 
will be measured as either t oi" i with equal probability) . 
Therefore, in the z direction, the spin intensity measured 
for each configuration in terms of these underlying states 
will be given by: 

^parallel , -^parallel , -^antiparallel 



J -^antiparallel -^antiparallel -^parallel /CC'\ 

^et,e; 7, \ \ — • \^^) 



Here /parallel and /antiparallel represent the total photoe- 
mission intensity given off by the states discussed in the 
main text; the single terms on the left represent the spin- 
polarized photoemission described in the main text (di- 
vided by two due to only half of the photons interacting 
with vni^ \ while the two terms on the right represent the 
non-spin-polarized photoemission (divided by four due to 
half photons interacting with mi^ , and half of those be- 
ing measured as either up or down). This dilution causes 
a decrease in the observed splitting compared to the in- 
trinsic one, which can be corrected by taking 



/parallel = (3/e;,et - l^et,ei)/2 = /t(1 - 2P^)/2 

/antiparallel = (3/et,ei " 1^01, et)/^ = ^t{\ + 2P®)/2 

(S9) 

to recover the intensities from each state; this assumes 
that all observed intensity is from only these two states 
with no background. Such a correction is shown in 
Fig. S3(f) and S4(f), as weh as Fig. 2(d) of the main text, 
and returns a larger splitting. 

Observing the uncorrected splitting at the different 
photon energies and energy ranges measured results in 
estimates of 83meV, 101 meV, and 110 meV for an aver- 
age of 98 ± 14 meV. However, this estimate must be low 
due to the experimental geometry and its effect as de- 
scribed above. Using the correction, as described above, 
on the same data sets results in estimates of 125 me V, 
184 meV, and 172 meV for an average of 160 ± 30meV. 
However, this estimate may be high, as it assumes that 
all intensity arises from these two states with no back- 
ground. These limits therefore represent upper and lower 
bounds of the intrinsic splitting between these two states, 
for an overall estimate of 130 ± 30meV. 




FIG. S5: (sfc • s-k ) between two states with opposite mo- 
mentum on the Fermi surface sheets in Sr2Ru04 for the {3 
(a,d,g), 7 (b,e,h), and a (c,f,i) bands at /c^ = (g-i), f (d-f), 
and TT (a-c). The location in the kx-ky plane is defined by 
the angle G for each band, as illustrated in the inset at lower 
right. 



Available spin-eigenstates around the Fermi surface 



Thus far we have examined only single-particle prop- 
erties, however, without making any assumptions about 
the pairing mechanism, we may also use our model to 
examine the spin-eigenstates available to a pair of elec- 
trons with net momentum - one electron at k and one 
at — /c. These eigenstates are shown in Fig. S5 by plotting 
the expectation value of s^-S-k around the Fermi surface 
sheets for all bands at three kz values. Pure, uncoupled, 
spins will form the familiar singlet /triplet eigensystem, 
which is visible at a single state with {s^ • S-k) = — f 
and a triplet of states with {si • S-k) = \] such a 
system can be seen to occur off the zone diagonal for 
the a band in Fig. S5(c,f,i). However, a singlet /triplet 
eigensystem is not available for either the /3 or 7 bands 
anywhere in momentum space, as seen in Fig. S5(a,d,g) 
and (b,e,h), respectively; instead the available eigensys- 
tem is observed to be either a doublet and two singlets, 
or two doublets, depending on the location in /c-space. 
In this regard, we also note that the introduction of a 
momentum-dependent pseudo-spin would not resolve this 
issue, since the unitary transformation a = U(k) a that 
in a momentum-dependent fashion links the spin and 
pseudo-spin for the one-particle eigenstates also enters 
in the spin operator, s = U^(k) s U(k), thereby leaving 
the expectation values unaltered. 
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